Clinical PDT dose dosimetry for pleural Photofrin-mediated photodynamic therapy

Abstract. Significance Photodynamic therapy (PDT) is an established cancer treatment utilizing light-activated photosensitizers (PS). Effective treatment hinges on the PDT dose-dependent on PS concentration and light fluence-delivered over time. We introduce an innovative eight-channel PDT dose dosimetry system capable of concurrently measuring light fluence and PS concentration during treatment. Aim We aim to develop and evaluate an eight-channel PDT dose dosimetry system for simultaneous measurement of light fluence and PS concentration. By addressing uncertainties due to tissue variations, the system enhances accurate PDT dosimetry for improved treatment outcomes. Approach The study positions eight isotropic detectors strategically within the pleural cavity before PDT. These detectors are linked to bifurcated fibers, distributing signals to both a photodiode and a spectrometer. Calibration techniques are applied to counter tissue-related variations and improve measurement accuracy. The fluorescence signal is normalized using the measured light fluence, compensating for variations in tissue properties. Measurements were taken in 78 sites in the pleural cavities of 20 patients. Results Observations reveal minimal Photofrin concentration variation during PDT at each site, juxtaposed with significant intra- and inter-patient heterogeneities. Across 78 treated sites in 20 patients, the average Photofrin concentration for all 78 sites is 4.98  μM, with a median concentration of 4.47  μM. The average PDT dose for all 78 sites is 493.17  μMJ/cm2, with a median dose of 442.79  μMJ/cm2. A significant variation in PDT doses is observed, with a maximum difference of 3.1 times among all sites within one patient and a maximum difference of 9.8 times across all patients. Conclusions The introduced eight-channel PDT dose dosimetry system serves as a valuable real-time monitoring tool for light fluence and PS concentration during PDT. Its ability to mitigate uncertainties arising from tissue properties enhances dosimetry accuracy, thus optimizing treatment outcomes and bolstering the effectiveness of PDT in cancer therapy.


Introduction
Photodynamic therapy (PDT) is an evolving treatment modality that has gained approval from the U.S. Food and Drug Administration for the treatment of various cancers, including microinvasive lung cancer, obstructing esophageal cancer, and non-small cell lung cancer.2][3][4] PDT offers valuable advantages, such as minimal invasion, low systemic toxicity, versatility, and repeatability.However, the optimal clinical outcomes of PDT have been hindered by the lack of a reliable dosimetry metric that accurately quantifies the effective treatment dose, which is essential for predicting and guiding PDT. 5 The challenges in achieving accurate dosimetry for PDT stem from its dynamic and complex nature, involving intricate interactions between light, the photosensitizer (PS), and tissue oxygen ( 3 O 2 ). 6In a typical PDT process, the PS is excited by specific wavelengths of treatment light, causing it to transition from its ground state to the excited singlet state.Subsequently, the PS undergoes intersystem crossing, leading to its transition to the triplet state.Type II PDT, the most clinically relevant process, occurs when the triplet state transfers energy to ground-state oxygen, 3 O 2 , generating singlet oxygen, 1 O 2 .[9][10][11] Conventional PDT dosimetry relies on a prescribed administered drug dose and the total light fluence, which represents the energy per unit area delivered by the end of treatment.3][14] Thus, the improvement of dosimetry methods is crucial for advancing the application of PDT.][17] Our research group has developed a dosimeter that enables the simultaneous measurement of light fluence rate and PS fluorescence during treatment.9][20] The measured PS fluorescence at the tissue surface allows for the calculation of PS concentration.Together with the measured light fluence, this information facilitates the determination of the PDT dose delivered to the superficial tissue of the pleural cavity.We have observed spatial heterogeneities in the delivered PDT dose on the cavity surface due to complex pharmacokinetics and tissue conditions. 14Consequently, it is crucial for the dosimetry system to cover the entire cavity surface rather than focusing on a single isolated point.To address this, the newly developed dosimetry system features eight dual-function channels, enabling the monitoring of eight different locations.The primary aim of this study is to provide updated information on PDT dose variations within the context of our pleural PDT clinical trial.Notably, the study expands its scope by increasing the number of patients from 8 to 20 and the inclusion of sites from 22 to 78. 14,21 Furthermore, our research involves a thorough and detailed examination of the obtained data.This includes a comprehensive reanalysis of fluorescence data to derive PS uptake, a thorough reevaluation of tissue optical properties through Diffuse Reflectance Spectroscopy (DRS) data, and a complete reassessment of optical property correction factors for determining PS uptake in both patients and phantoms.Through these efforts, we address previously existing gaps in the literature, providing a more comprehensive and updated understanding of the subject matter.In addition, we aim to incorporate an optical infrared (IR) navigation system to provide comprehensive dosimetry guidance throughout the pleural PDT procedure.Although existing real-time dosimetry methods primarily focus on in-vivo light fluence measurements, the use of innovative instruments makes real-time PDT dose evaluation feasible.This study provides a summary of the current state of the art in multichannel PDT dose measurements for clinical applications and explores the various potentials of this versatile dosimetry system.

Clinical Trial and Concurrent Measurements
Patients with pleural mesothelioma were enrolled in a phase II/III randomized clinical trial at the Hospital of the University of Pennsylvania after providing informed consent.A subset of patients underwent lung-sparing surgical resection combined with interoperative Photofrin-mediated PDT, and the remaining patients received surgery alone.Further details regarding the PDT technique and treatment protocol can be found elsewhere. 22For patients receiving PDT, Photofrin (Pinnacle Biologics, Chicago, Illinois) was administered intravenously at a dose of 2 mg per kg of body weight, 24 h prior to PDT.PDT was performed immediately following surgery using a KTP-pumped dye laser (model 630 XP, Laserscope, Inc., San Jose, California) with 632 nm light.The prescribed total light fluence was set at 60 J∕cm 2 .The instantaneous light fluence rate was measured using eight isotropic detectors (Medlight SA, Ecublens, Switzerland) sutured at the apex, posterior chest wall (PCW), anterior chest wall (ACW), posterior sulcus (PS), anterior sulcus (AS), posterior mediastinum (PM), pericardium (Peri), and diaphragm (Diaph) 23 on the pleural cavity surface (see Fig. 1).The cumulative light fluence was simultaneously calculated until the prescribed light dose was achieved at all locations.In addition to the light fluence rate, the Photofrin fluorescence excited by the treatment light was monitored using the same eight isotropic detectors.The light was delivered via a bare optical fiber embedded in a modified endotracheal tube.To facilitate light scattering, the treatment tube and pleural cavity were filled with 0.1% Intralipid.The physician sequentially moved the light source along the inner surface of the cavity during PDT until the prescribed light dose was reached. 24ccurately quantifying fluorescence emission in vivo is challenging due to the varying optical properties of the surrounding tissue.To mitigate this effect, diffuse reflectance measurements were obtained on the inner surface of the pleural cavity before and after PDT treatment for optical property correction.A custom-made fiber optic-based contact probe (FiberOptic Systems, Inc., Simi Valley, California) was employed, consisting of one source fiber connected to an air-cooled quartz-tungsten-halogen lamp (Avalight HAL-S, Avantes, Inc., Louisville, Colorado) and nine detection fibers spaced at distances ranging from 1.4 to 8.7 mm from the source fiber.The detection fibers captured the diffuse reflected light, which was then directed to a spectrograph and recorded on a charge-coupled device (CCD)-based camera sensor (InSpectrum 150, Roper Scientific, Princeton, New Jersey) for radially resolved diffuse reflectance measurement.Background signals were concurrently measured for subsequent subtraction.The acquired spectra were fitted using a nonlinear fitting algorithm implemented in the MATLAB programming environment (The MathWorks, Natick, Massachusetts) to extract the tissue optical properties at each measured location.Further information regarding the probe design and fitting algorithm can be found elsewhere. 21,25g. 1 Point cloud data of a reconstructed pleural cavity, with isotropic locations marked, were generated using our scanning system.

PDT Dose Dosimeter Instrumentation
The comprehensive multichannel dosimetry system incorporated the use of eight isotropic detectors that were sutured onto the interior pleural surface.These detectors were connected to individual channels within the dosimetry system.Figure 2 provides a front view of the multichannel dosimetry system and a schematic diagram of its internal setup.The system enclosure contained a total of 16 channels [Fig.2(a)].The internal structure of the equipment is shown in Fig. 2(b).The bottom row was comprised of channels 1 to 8, with each channel split into two fibers through 1-to-2 bifurcated fibers (Ocean Optics, Dunedin, Florida) internally inside the box.The bifurcated fibers of each channel were connected to the photodiodes and one spectrometer (Exemplar, B&W Tek, Inc., Newark, Delaware) for the simultaneous monitoring of the fluence rate of the treatment light and Photofrin fluorescence.To block the treatment light, a long-pass filter was installed for each channel in conjunction with the spectrometer.Due to space constraints, the remaining channels (channels 9 to 16) do not have the bifurcations, so they cannot measure the PS uptake and were solely connected to the dosimetry system.The PDT dose dosimetry instrument always had 16 channels for light dosimetry and was incrementally expanded to 8 channels to measure the PS uptake and light fluence rate simultaneously.
The spectrometers were utilized for measuring the Photofrin fluorescence, and the photodiodes served as monitors for the fluence rate of the treatment light.Recently, the dual-function channels were expanded from four to eight, allowing for a more comprehensive measurement of the entire cavity.The spectrometer had a wavelength range of 200 to 1050 nm and a resolution of 0.42 nm using a diode array with dimensions of 2048 × 1 elements and an element size of 14 μm × 200 μm.The achieved spectral resolution was 0.47 nm.The isotropic detector directly measured the light fluence rate at the surface of the cavity.However, it should be noted that this measurement might differ from the intra-tissue light fluence rate due to variations in tissue optical properties.To ensure the accurate measurement of the transmitted fluorescence (signal after 633 nm), long pass filters (Semrock, Inc., Rochester, New York) were employed to block the treatment light.No filtration was required for the treatment-light signal prior to reaching the photodiodes. 26In this study, all eight isotropic detectors were utilized and connected to channels 1 to 8 to simultaneously measure the fluence rate of the treatment light and Photofrin fluorescence throughout the pleural cavity.

Spectroscopy
Fluorescence spectra were collected for eight strategic locations within the pleural cavity wall, as described in Sec.2.1 and shown in Fig. 1, using eight single-channel CCD spectrometers.To evaluate the concentration of Photofrin, the raw fluorescence spectra obtained during PDT were corrected for the spectral response of each spectrometer.Subsequently, a single value decomposition (SVD) fitting algorithm was applied to analyze the spectra using a basis spectrum [Fig.3(a)]. 27The basis spectrum consisted of a laser component and a Photofrin fluorescence component, which were established through extensive phantom studies.The laser component served as a reference for the excitation light intensity.Background spectra were measured and subtracted.In the SVD algorithm, a 21-term Fourier series was incorporated to account for any unknown spectroscopic components, such as ambient room light.However, its weight in the fitting routine was lower than that of the measured fluorescence emission components.Throughout the treatment, each channel captured hundreds of fluorescence spectra, which were individually fitted using the SVD algorithm.The peaks observed around 675 nm corresponded to the fluorescence signal of Photofrin detected by the isotropic detectors.Slight shifts in the peak wavelengths were observed across different channels, attributed to the distinct long-pass filter characteristics of each channel.To accommodate this, during SVD fitting, adjustments were made to the basis spectrum to align with the shifted peaks for different channels.By fitting the raw measured spectra with the basic components, unitless SVD amplitudes were obtained.The laser components enabled the elimination of variations in the excitation light, allowing for the quantification of the local Photofrin concentration.A linear relationship was found between the SVD amplitudes and the local Photofrin concentration, enabling the quantification by multiplying with a constant.The determination of this constant was demonstrated in a previous study. 26The SVD fitting method utilized a large number of fluorescence spectra measured throughout the treatment at each location, leading to a significant reduction in uncertainties.The SVD fitting results were later adjusted to accommodate distortions arising from variations in optical properties, as elaborated in detail in Sec.2.4.

Optical Property Correction
][30][31] To account for the influence of different optical properties on the detected fluorescence signal, a correction factor based on absorption coefficient (μ a ) and scattering coefficient (μ 0 s ) was introduced. 24,32This correction factor, CF, was determined using a series of tissue-simulating phantoms with varying optical properties (μ a ¼ 0.06 to 0.7 cm −1 and μ 0 s ¼ 4 to 20 cm −1 ) while keeping the Photofrin concentration constant at 3 mg kg −1 .Intralipid (Fresenius Kabi, Germany) was used as a light scatterer, and ink (Parker ® Quink ® ) served as the light absorber in the phantoms.The values of μ a and μ 0 s represent the absorption and reduced scattering coefficients, respectively, at the emission wavelength of 630 nm, corresponding to the Photofrin fluorescence.The fluorescence signals from these phantoms were measured using the 8-channel PDT dosimeter.The raw fluorescence spectra obtained were then fitted using the SVD fitting algorithm described in the previous section, yielding an SVD amplitude (A SVD ) for each phantom, as shown in Fig. 3(c).An empirical formula fitting the experimental data was employed to calculate the optical property correction factor, CF, from a given tissue optical properties ðμ a ; μ 0 s Þ to that of the reference tissue optical properties (μ a;ref ¼ 0.32 cm −1 and μ 0 s;ref ¼ 10 cm −1 ): where ; t e m p : i n t r a l i n k -; e 0 0 2 ; 1 1 4 ; 2 4 4 Monte Carlo calculations were conducted to determine the fitting function and its parameters. 26,33A different set of tissue-simulating phantoms, featuring varied Photofrin concentrations ranging from 0 to 5 mg∕kg, was employed to create a calibration curve for the A cor and the absolute Photofrin concentrations.The resulting calibration curve is shown in Fig. 3(d).This calibration curve enables the determination of absolute Photofrin concentrations by utilizing the optical property correction factor (A cor ).

Real-Time Scanning and Navigation System
To enhance the prediction and efficacy of PDT, our research group has developed a comprehensive real-time scanning and navigation system, in addition to the light dosimetry system.This system consists of a commercial optical infrared (IR) camera (Polaris, NDI, Waterloo, Canada), trackable wands, and a 3D scanner, as shown in Fig. 4. The hand-held 3D scanner (Structure Core, Occipital, Boulder, Colorado) is utilized for rapid and precise capture and reconstruction of the pleural cavity, enabling the identification of the target surface for real-time calculation of light fluence distribution during PDT. 34The IR navigation system comprises stereo cameras, a modulated laser source (wavelength 850 nm), and a treatment wand equipped with the PDT light source and nine passive markers.The cameras track the position of the light source in real time by detecting the light reflected from the markers on the wand at a rate of 20 to 60 Hz.This navigation system has been successfully employed in clinical trials for both HPPH-and Photofrin-mediated pleural PDT.The collected data from these trials have been utilized for post-treatment analysis, providing valuable insights into the treatment process and facilitating outcome prediction. 35By incorporating this innovative scanning and navigation system alongside the dosimetry system, we can provide more comprehensive real-time guidance for PDT, leading to improved treatment efficiency.
To provide real-time feedback on the light fluence distribution across the entire inner surface of the cavity, a newly developed graphical user interface (GUI) was implemented. 36,37This GUI enables the visualization of the light fluence distribution in real time.Simultaneously, the light dosimetry system was utilized to measure the light fluence for verification purposes.It should be noted that the light fluence at each point of interest on the inner surface of the cavity is the sum of both the primary component and the scattered component of the treatment light. 22The GUI provides a comprehensive and real-time representation of the light distribution, allowing for enhanced monitoring and optimization of the PDT procedure.The primary component of the light fluence rate (φ) is calculated as E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 0 3 ; 1 1 7 ; 2 1 1 ϕ ¼ S 4πr 2 ; ( where S is the source power (mW) and r is the distance (mm) between the point light source and the cavity surface.For more accurate calculation, scatter light fluence rate is included; it is calculated as E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 0 4 ; 1 1 7 ; 1 4 1 where b is an empirical constant. 23,37The ultimate goal is to integrate this innovative system with the newly developed 8-channel dosimetry system, enabling the provision of real-time PDT dose information for the pleural cavity during PDT.
Fig. 4 Instruments consist of the newly developed real-time navigation system.The optical infrared (IR) navigation system is made up of (a) an IR camera and (b) a treatment wand with a laser fiber inserted.The novel scanning system is made up of (c) a 3D scanning device for capturing the pleural surface.
3 Results and Discussion

Tissue Optical Properties and Correction Factors
In this study, in vivo diffuse reflectance measurements were conducted at 52 sites within the pleural cavities of 11 patients.Initially, the PDT dose dosimetry system had four channels capable of simultaneously measuring light fluence rate and fluorescence.However, for the latest two patients, the system was expanded to eight channels, covering all eight predetermined locations within the pleural cavity for each patient.In this paper, the presented data for the four-channel system differ from previously reported results obtained using the same system.All raw data underwent thorough validation and were re-analyzed specifically for this study.Table 1 summarizes the tissue absorption and reduced scattering coefficients (μ a and μ 0 s ) at the excitation wavelength of 630 nm for all measurement sites where diffuse reflectance spectra were obtained, along with the corresponding correction factors.Significant heterogeneity in both absorption and reduced scattering coefficients can be observed both within and among patients.
The maximum variation of μ a and μ 0 s is 8 times and 2.8 times, respectively, among all sites with a mean (SD) of μ a ¼ 0.31 (0.13) and μ 0 s ¼ 9.9 (2.4) cm −1 , respectively.The variations in tissue optical properties can be attributed to factors such as thermal damage to the surrounding tissue induced by electrosurgery during tumor resection and the diverse characteristics of cellular content and chromophores throughout the tissues. 38These factors are uncontrollable during the treatments, necessitating post-treatment correction.To provide a comprehensive understanding and address the influence of optical properties, correction factors derived from measurements taken at 78 sites across 20 patients at the treatment wavelength are illustrated in Fig. 5.This includes data from earlier cases using 2-or 4-channel PDT dose dosimetry systems.CFs were computed for individual sites using Eq. ( 1) to adjust the measured fluorescence for distortions arising from variable optical properties at different sites and among different patients.This correction is essential for ensuring the accurate quantification of PS concentration.The CFs range from 0.54 to 2.57 across the 78 sites.The mean CF value is 1.0 (indicated by the dashed line), and the median CF value is 1.08.The values for each site are presented in Table 1.The maximum variation among all sites is 4.7 times, and the maximum variation among sites within individual patients is 1.8 times.

Temporal and Spatial Distribution of Photofrin and PDT Dose
Figure 6 shows the temporal changes in local Photofrin concentrations for 16 sites in the most recent two cases, in which all eight dual-function channels were utilized.The discrete data points in the graph were obtained from measured fluorescence spectra using the SVD fitting method described earlier.It is worth noting that some data points exhibit larger variations, which can be attributed to the fluctuations in the incident treatment light fluence rate (ranging from 0 to 1200 mW∕cm 2 ).This variation stems from the movement of the light source during PDT, which introduces errors in the measurement of Photofrin fluorescence excited by the treatment light.
To facilitate a clearer observation of the trend of Photofrin variation during the treatment, the fitted data points for each channel were smoothed by calculating the average of all data points every 10 min of treatment time.The smoothed results are depicted by the solid curves in the graph.The relatively small temporal fluctuations in Photofrin concentration throughout the treatment for each site suggest that, despite the significant variations observed between different locations, the changes in Photofrin uptake over time within each site remain relatively stable.This phenomenon suggests that the uptake of Photofrin is not significantly influenced by the treatment itself but exhibits variation across different locations.The observed differences in optical properties (μ a and μ 0 s , tabulated in Table 1 for each patient) at these locations can be attributed to the inherent complexity of the tissue environment.This observation aligns with the substantial variation in CFs (see Fig. 5), calculated based on the tissue optical properties at the sites, as described in Sec.3.1.Importantly, there is no significant photobleaching of Photofrin during PDT, corroborating findings from previous studies. 26he mean Photofrin concentrations (shown on the left axis) and PDT doses (shown on the right axis) obtained from all 78 sites in the pleural cavities of 20 patients are shown in Fig. 7.The average (SD) and median Photofrin concentrations for all 78 sites are 4.98 (2.38) μM and 4.47 μM, respectively.The PDT dose delivered to each site is defined as the product of the local Table 1 Summary of absorption coefficients, scattering coefficients, and correction factors for 20 patients.The light fluence on the surface is the same at 60 Jcm −2 for all patients.Information about patients 07, 08, 12, 14, 16, 17, 18, and 20 can be found elsewhere and were reprocessed in this study. 26

Patient
Site* Optical properties  E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 0 5 ; 1 1 4 ; 2 5 0 PDT dose μMJ The differences in effective PDT doses delivered to different sites were influenced by the heterogeneities in Photofrin uptakes along the pleural cavities.The delivered PDT doses for each site are shown on the same plot (Fig. 6) with a secondary vertical axis on the right to indicate the corresponding values.The average (SD) and median PDT doses for all 78 sites are 493 (236) μMJ∕cm 2 and 442.79 μMJ∕cm 2 , respectively.The maximum variation among all sites is 9.8 times, indicating significant differences in the delivered PDT doses across the pleural cavities.Similarly, the maximum variation among sites within each patient is 3.1 times, highlighting the variability in PDT doses within individual patients.
The substantial variations observed within and among patients align with previous studies and underscore the importance of considering both Photofrin concentration and PDT dose in the development of an effective dosimetry system.Traditional approaches that solely rely on monitoring light fluence may not provide comprehensive information to ensure uniform PDT doses across all treatment sites.One potential limitation in real-time monitoring of Photofrin concentration is the time-consuming post-treatment fitting and analysis of measured fluorescence spectra.However, this study demonstrates relatively stable Photofrin concentrations throughout the time course of PDT, enabling the utilization of data acquired at the beginning of treatment for comprehensive real-time guidance throughout the rest of the procedure.Therefore, the innovative 8-channel PDT dose dosimeter holds the potential to assist physicians in delivering more uniform PDT doses to all sites, ultimately improving the outcomes of pleural PDT treatment.

Real-Time Navigation System
The preliminary results for light delivery efficiency with and without real-time guidance using the scanning and navigation system are depicted in Fig. 8.A pleural phantom (see Fig. 1), created using a 3D printing technique and based on a CT scan of an actual human pleural cavity, was utilized in our real-time navigation study.To simulate clinical conditions more accurately, six isotropic detectors were placed at predetermined locations on the inner surface of the pleural phantom (Apex, PCW, ACW, PS, PM, Peri) and connected to the light dosimetry system.The positions "Diaph" and "AS" were not included due to the positioning of the opening on the phantom, which allowed the scanner and treatment light source access.By comparing the results, it can be observed that employing the scanning and navigation system leads to a more uniform overall light distribution compared with PDT treatment using detectors alone.Although the prescribed light doses are accurately achieved at the crucial positions where the isotropic detectors are fixed, there are some remaining "cold spots" with lower light fluence.The innovative real-time navigation system effectively addresses this issue and enables the achievement of uniform light fluence distribution by the end of the treatment.The infrared (IR) navigation system has been utilized in previous clinical pleural PDT trials to track the motion of the light source during PDT and evaluate the uniformity of light fluence distribution across the entire pleural cavity after treatment.This serves as a baseline for the application of the updated navigation system, which incorporates rapid scanning and real-time feedback capabilities.The integration of these advancements can further enhance the accuracy and effectiveness of pleural PDT procedures. 39he real-time light fluence calculation of the novel system was validated by comparing the calculation with the data measured by all six detectors used in the phantom studies with different optical properties.To account for the general scattering inside the cavity and reduce the percentage error, a constant "b" was introduced based on previous studies.In this study, the scatter component was applied after treatment.The percentage error from the measured light dose at the end of treatment with real-time calculated light fluence using the primary component only [Eq. ( 1)] and the primary and scattering component together [Eq.( 2)], along with the respective "b" values used, are summarized in Table 2.The variation in "b" values was found to be small among previous clinical cases, and this study confirms the same observation.After applying the scatter component, the range of the percentage error in light fluences is between 0.9% and 12.8%, which is comparable to the findings of previous studies (ranging from 0.7% to 15.4%). 23his study demonstrates the improvement in light delivery efficiency achieved by utilizing the real-time navigation system, and the algorithm employed is validated.The validation process affirms the accuracy and reliability of the real-time light fluence calculation, establishing its suitability for application in pleural PDT procedures.This establishes a foundation for the development of a comprehensive real-time guidance system, integrating the 8-channel dosimetry system with the navigation system.Through the combination of real-time calculated light fluence data for the pleural cavity and simultaneously measured PS concentrations at eight crucial sites, our objective is to furnish physicians with real-time PDT dose information.This information serves as a guidance tool for optimizing light delivery during PDT, with the overarching goal of further advancing treatment efficiency.In this study, a real-time 8-channel PDT dose dosimeter was developed and utilized for Photofrin-mediated pleural PDT.The mean PDT dose among the 20 patients included in the study was found to be 493.17μMJ∕cm 2 .However, it was observed that, with the same administered Photofrin dose and total light dose, PDT doses could vary significantly, ranging up to 980% among different patients and 310% within the same patient.This variability implies that, relying solely on light fluence as the treatment guidance during PDT leads to a non-uniform distribution of PDT dose throughout the pleural cavity, potentially impacting the overall treatment outcome.Recognizing PDT dose as a crucial dosimetry parameter, it can serve as a superior guiding tool for physicians during PDT.To address this issue and improve the consistency of PDT doses, the 8-channel dose dosimeter was developed, providing a means to initiate PDT dose-mediated pleural PDT.The system incorporated commercial tracking and scanning tools along with customized software to enable real-time feedback and monitoring.The scanning and navigation system efficiently identified the target surface and facilitated quick registration for real-time light fluence distribution calculation during PDT.The developed system has been successfully applied in phantom studies, demonstrating its ability to provide reliable real-time feedback on the light source position and 2D light fluence distribution.
The significance of the 8-channel system lies in its potential to provide informative real-time guidance at all sites currently used in the clinical trial.By addressing variations in PDT doses and ensuring uniform dose distribution, the system offers a promising avenue for starting a new clinical trial utilizing PDT dose as dosimetry metrics.Future development will prioritize the integration of the 8-channel dosimetry system with the scanning and navigation system to guide light delivery during PDT.This integration aims to guarantee the consistent distribution of PDT dose, thereby further elevating the effectiveness of cancer treatment outcomes.

Disclosures
Theresa M. Busch received support from Simphotek and personal fees from Lumeda and IBA outside the submitted work.Keith A. Cengel received support from Simphotek.All other authors declare no commercial conflicts of interest related to this article.

Fig. 2
Fig. 2 (a) The front view of the 8-channel PDT dose dosimetry instrument and (b) the schematic diagram of the system setup.The bifurcated fibers are connected to channels 1-8 internally, enabling the simultaneous measurement of light fluence and PS uptake.A long-pass filter is employed for each channel, leading to a slight variance in the peak wavelength of the measured spectrum.To ensure independent functionality, distinct long-pass filters are employed for each channel.

Fig. 3
Fig. 3 (a) The basis spectrum consists of a laser component, a Photofrin fluorescence component, and a Fourier component; (b) example of measured raw fluorescence spectra from 8 channels of patient #53.The peaks at 675 nm arise from fluorescence of the isotropic detector.Slight shifts in the peak wavelength occur as a result of using different long-pass filters for each channel.(c) Fluorescence SVD amplitudes for Photofrin in tissue-simulating phantom experiments with different optical properties (μ a ¼ 0.06 − 0.7 cm −1 and μ 0 s ¼ 4 to 20 cm −1) with a constant Photofrin concentration of 3 mg kg −1 , measured in phantoms using 8-channel PDT dosimeter (solid lines) and Aðμ a ; μ 0 s Þ fits using Eq.(2) (dashed lines).The correction factor (CF) and the corrected Photofrin SVD amplitudes were obtained using Eqs.(1) and (2).(d) Photofrin concentration calibration curve.
Photofrin concentration ([PS], in μM) and the delivered light dose.The treatment for each site was halted when the measured total light fluence of 60 J∕cm 2 was reached, ensuring a consistent total fluence across all subjects:

Fig. 5
Fig. 5 CF at 78 different sites in the pleural cavities of 20 patients.(a) Patient #07 to #29.(b) Patient #32 to #52.The mean value is represented by the dashed line.

Fig. 6
Fig. 6 Temporal changes of Photofrin concentrations measured from 16 sites in the pleural cavity of the 2 most recent patients (#052 and #053) during the PDT treatments.To convert from mg kg −1 to μM, 1 mg kg −1 ¼ 1.65 μM can be used.

Fig. 7
Fig. 7 Mean Photofrin concentrations and PDT dose delivered to 78 different sites in the pleural cavities of 20 patients.(a) Patient #07 to #29.(b) Patient #32 to #52.The mean value is represented by the dashed line.

Fig. 8 Table 2
Fig. 8 Comparison of light fluence map at the end of treatment with six detector positions labeled.(a) Without real-time fluence map guidance.(b) With real-time fluence map guidance.The overall light fluence distribution presented for the whole surface area was calculated based on the navigation system and Eqs.(3) and (4).Light fluence data were measured by the isotropic detectors at the locations marked with "X" to verify the calculations.